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The gluon distribution of a nucleon propagating through a nucleus can change 
depending on how far the penetration depth is. It is a nonperturbative process 
that we describe by an evolution equation. The kernel of the integral equation is 
to be determined by a phenomenological study of J/ip suppression in pA collision, 
treated as a perturbative process. The data of E866 on a(xp) shows significant 
dependence on xp at large xp. It presents a feature that has not been explained 
by any dynamical model. We show that gluon depletion is a simple mechanism 
that can account for it. The result has far-reaching implications on the role of 
partons in nuclear collisions. 



The conventional reason for studying J/ip suppression in heavy-ion colli- 
sions is that it might reveal evidences for the existence of quark-gluon plasma 
which enhances the deconfinement of the charmonium states. Thus the the- 
oretical attention is usually focused on the effect of the medium on the cc 
states that are created by gluons. Most experimental work that reveal the 
suppression phenomenon in nucleus-nucleus (AB) collisions have been done at 
CERN-SPS, where the measurement is limited to the central rapiditji.jjcgion. 
Indeed, impressive results have stimulated a great deal of excitement aati. 

The emphases in this paper are unconventional □. First, we consider pA 
collisions for which we know that quark-gluon plasma is not likely to form. 
Second, we examine the J/ip at all positive xp, not just around xp ~ 0. 
Third, we focus on the possibility that the gluon distribution is significantly 
modified as the projectile traverses the target nucleus. Note that we use the 
word "projectile" instead of "proton" because what goes the nucleus can be so 
distorted from the incident proton that a realistic description of the "projectile" 
may well be just a flux of partons. Indeed, the significance of this work is in 
sowing the seeds of suspisicion that the conventional notions of wound nucleon 
and participating nucleons may be outmoded. 

The stimulus for this line of investigation is the FNAL-E866 data on J/ip 
suppression in pA collisions at large xp The cross section for J / ip production 
is found to depend on A and if as 

v p a->j/4,(xf) = o- pp ^j/^{x F )A a{xF) . (1) 

In Ref. 5 the values of a (xp) are tabulated; moreover, a fit of a {%f) by an 
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analytic formula is given: 



a(x F ) = 0.96 (1-0.0519^-0.338:4) . (2) 

Evidently, there is a sizable amount of suppression at large xp. The issue is 
how this enhanced suppression can be understood theoretically. 

The conventional approach, as stated in the beginning, is to study the 
effects of the medium on the produced. cp state. That has been carried out by 
at least three groups of investigatorsOlIm. The consensus is that if the cc state 
is to break up quickly at high xp , the time involved is too short (~ 0.02/m/c) 
to be feasible in the usual nuclear medium. The conclusion is that there is 
no mechanism in hadronic absorption of the final state capable of explaining 
the Xf dependence. If we use H (A) to denote that hadronic absorption effect, 
it is reasonably safe at this stage to assume that it isJnsensitive to xp. The 
canonical description of H(A) is the exponential formEl 

H(A) = exp [-paz(A)] , (3) 

where p is the average nuclear density, a the absorption cross section, and z(A) 
is the average path length in A that the cc state traverses. 

If the final state dynamics cannot account for the xp dependence, then 
the only alternative is in the initial state. The subprocess oi g + g — >c + c is 
cancelled in the ratio R = <j p a^,j /ipl ' Aa pp ^j There are two types of effects 
on the initial state, both having to do with the gluon flux that leads to the 
subprocess g (xi) + g {x%) — > cc(xp). Here g (x\) is the gluon with momentum 
fraction x\ of the initial proton, g [x^] is the gluon with momentum fraction 
X2 of some nucleon in the target nucleus, and xp = x\ — X2- The two effects 
are gluon depletion associated with g(xi), and nuclear shadowing associated 
with g (22). Denoting the two by G [x\, A) and (x2, A), respectively, we can 
write the ratio R, defined earlier, by 

R(xf,A) = G(xi(x f ),A)N(x 2 (xf),A)H(A) . (4) 

From ([!]) we know empirically 

R(x F ,A) = A a{XF) - 1 . (5) 

Since there is some independent information on N (X2 ,A), to be discussed 
shortly, it is clear from (||)-((|) that the behavior of G{x\,A) is highly con- 
strained by what is already known phenomenologically. Our goal is to de- 
termine G{x\, A) on the basis of a sensible evolution equation for the gluon 
distribution. 
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On nuclear shadowing there exists now extensiye quantification of the 
quark and gluon distribution functions in nuclei EjO. They are determined 
by analyzing the deep inelastic scattering and dilepton production data of nu- 
clear targets at high Q 2 on the basis of DGLAP evolutionlij. IDie results are 
given in terms of numerical parameterizations (called EKS98 ELl) of the ratio 
Nf (x,Q 2 ) = fi/A (x,Q 2 ) I fi (x,Q 2 ), where /j is the parton distribution of 
flavor i in the free proton and fi/A is that in a proton in a nucleus A. we shall 
be interested in the ratio for the gluon distribution only at Q 2 = lOGeV 2 , 
which corresponds to producing cc near the threshold. We denote that ratio of 
the gluon distributions by N(x2,A), which is what appears in (Q). A simple 
formula was found that can provide a good fit of the EKS98 results; the details 
of which will appear in Ref.13. Here, we simply state the formula, used already 
in Ref.4. 

N(x 2 ,A) =A^ , (6) 

where 

0(£(x 2 )) = £ (0.0284 + 0.0008£ - 0.0041£ 2 ) , (7) 

and 

£ = 3.912 + £nx 2 . (8) 

For £ > 0, corresponding to x-2 > 0.02, (3 is positive and the region is usually 
referred to as anti-shadowing. The region of X2 that is relevant to our study 
of J /ip suppression at all positive values of xf turns out to straddle x 2 = 0.02 
and thus includes both shadowing and anti-shadowing. Using (^|) and (|^) in 
(Q) we have 

G(xi {x F ),A)H{A) = 4oOw)-0(«to(*F))-i , (9) 

whose RHS is now regarded as known. 

G{x\,A) is the ratio of the gluon distribution of a proton penetrating 
a nucleus, g(x\,A), to that in a free proton, g(xi,0). To be more precise, 
g (xx, A) should be labeled as g (xi,z(A)), where z{A) is the average distance 
that a proton propagates in the nucleus A before its gluon at x\ interacts with 
a gluon in the target at X2 to produce the cc pair. The modified distribution 
g(xx,z(A)) is, of course, independent of what hard subprocess takes place at 
x%; the label z(A) merely denotes how far the gluons have penetrated the 
nucleus A when that subprocess does occur. Let us then write explicitly 

G(x 1 ,A)=g(x 1 ,A)/g(x 1 ,0) . (10) 
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Using this in it is convenient to consider another function J (x\, A), defined 
by 



J(x 1 ,A)=g(x 1 ,A)H(A)=g(x 1 ,0)A'*-P- 1 . (11) 

For the gluon distribution in a free proton, it is sufficient for us to use the 
canonical form 

S(a:i,0) = ; g (l-xi) 5 , (12) 

since it is only the deviation from that form that is of interest; even the co- 
efficient go is unimportant for it will be cancelled out in the ratio (|o|). The 
significance of J (xi, A) is that it is completely known by virtue of the RHS of 
(pi]). Thus, in principle g(xi,A) is known phenomenologically. However, we 
need some theoretical input to give g (xi, A) a suitable analytical form. 

Since the effect of a nuclear target on the projectile gluon distribution is 
highly nonperturbative, there is no reliable analytical way to treat the problem 
of gluon depletion from first principles. Nevertheless, we propose an evolution 
equation in the spirit of DGLAPL3 except that we replace £nQ 2 by penetration 
length z in a nucleus. For the change of g(x, z), as the gluon traverses a distance 
dz, we write 

d f 1 dx' , . . / x \ 

d-z g(x > z) = J x — ' (13) 

where Q(x/x') describes the gain and loss of gluons in dz, but unlike the 
splitting function in pQCD, it cannot be calculated in perturbation theory. 
Eq. ( p~3| ) represents an approximation that ignores the quark channel, which 
should be included in a more complete treatment. Since Q {x/x') is unknown, 
we determine it phenomenologically from our knowledge of J(x,A) through 
0- 

Eq. (|13j) can easily be solved if it is put in the form of its moments. Define 
g n (z)=f dxx n - 2 g(x,z) (14) 

J X 

and similarly Q n . Then from the convolution theorem follows 

dg n (z)/dz = g n {z)Q n , (15) 

which yields 

9n(*) =9n(0)e zQ " . (16) 
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If we also take the moments of the left half of (pi]), we have J n (z) — g n (z)H(z), 
since H(A) is independent of x%. Using (|3|) and (|l6|), we obtain 

K n (z)=£n[J n (z)/g n (Q)]=z(Q n -pa) . (17) 

Since the quantity in the middle above is known, we can determine Q n . It 
turns that there are some subtleties that render the determination less than 
straightforward. For details the reader is referred to Refs. 4 and 13. 

Once Q n is known, we can use ( fL6| ) to calculate g n (z)> which in turn de- 
termines g{x\,z). The result can best be exhibited by G(xi,z). Using z to 
denote z(A), the average path length in A for the cc state to be produced, we 
show the behaviors of G(x\,A) for two representative values of A in Fig. 1. 
Evidently, there is an appreciable degree of suppression of the gluon distribu- 
tion function at high xi, and a small enhancement at low X\. The implication 
is that a gluon with high momentum can, upon interaction with the nuclear 
target, split up into gluons with lower momenta. 

The significance of our finding about the modification of the gluon distri- 
bution goes beyond the J /if) suppression problem itself, since it would revise 
the conventional thinking concerning the role of partons in nuclear collisions. 
The usual procedure for calculating hard processes is to use the parton distri- 
bution in the free proton as input for the hard subprocess. That procedure 
is based on factorization, which is obviously invalid if significant depletion of 
gluons takes placed in the nucleus. Furthermore, the redistribution of gluons 
as they propagate through the nucleus casts doubt on the plausibility of the 
notion that a penetrating nucleon can be identified as a recognizable entity, 
even if wounded. If so, the validity of the wounded-nucleon model that is con- 
ventionally used in nuclear collisions is now subject to serious reexamination. 
Perhaps even the meaningfulness of the discrete counting of the number of 
participants is questionable. All these issues need further investigation in light 
of our present finding. 

A way to confirm gluon depletion is to measure the suppression of strangeness 
production, either open or closed ss states, at large x p . The general strangeness 
enhancement can only partially be affected by the small gluon enhancement 
at small x%, as can be seen in Fig. 1. That small enhancement does not lead 
to any significant increase in the cc state either, so that the known J/ip sup- 
pression in the central region of AB collisions should be unaffected. A rather 
intriguing question is the behavior of J/ip suppression in pA collisions in the 
xp < region. Experimentally, that can be investigated in the inverse collision 
processes of nuclei on fixed proton target, or in colliding beams as at RHIC. 
Apart from the characteristics of parton redistribution in nuclei that will play 
a strong role in this problem, a new feature of gluon depletion can possibly 
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arise and be tested. In the proton rest from one has a row of nucleons moving 
the same direction, but the gluons belonging to the rear part of the row may 
be depleted due to their interactions with the slow gluons in the front part of 
row that are set free by their interactions first with the target proton. This 
type of effect, called nonlinear depletion, was first studied in Ref. 14, but could 
not be independently tested. But in the xp < region of pA collisions, that 
would be be a unique opportunity to study such effects. 
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Figure Captions 



Fig. 1 The ratio G(xi,A) of gluon distributions showing the effects of gluon 
depiction. 
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